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Abstract:  In  the  present  work,  we  use  the  homotopy  perturbation  method  (HPM)  to  solve  the  Newell- 
Whitehead-Segel  non-linear  differential  equations.  Four  case  study  problems  of  Newell- Whitehead- 
Segel  are  solved  by  the  HPM  and  the  exact  solutions  are  obtained.  The  trend  of  the  rapid  convergence 
of  the  sequences  constmcted  by  the  method  toward  the  exact  solution  is  shown  numerically.  As  a result 
the  rapid  convergence  towards  the  exact  solutions  of  HPM  indicates  that,  using  the  HPM  to  solve  the 
Newell-Whitehead-Segel  non-linear  differential  equations,  a reasonable  less  amount  of  computational 
work  with  acceptable  accuracy  may  be  sufficient.  Moreover  the  application  of  the  HPM  proves  that  the 
method  is  an  effective  and  simple  tool  for  solving  the  Newell-Whitehead-Segel  non-linear  differential 
equations. 
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INTRODUCTION 

Recently  lots  of  attentions  are  devoted  toward  the  semi-analytical  solution  of  real-life  mathematical 
modeling  that  is  inherently  nonlinear  differential  equations.  Most  of  the  nonlinear  differential  equations  do  not 
have  an  analytical  solution.  The  idea  of  the  homotopy  perturbation  method  was  first  pioneered  by  (He,  1999) 
Later  the  homotopy  perturbation  method  (HPM)  which  is  a semi-analytical  method  is  applied  to  solve  the  non- 
linear non-homogeneous  partial  differential  equations  (He,  2005;  Yildirim,  2009;  He,  2006;  He,  2006;  He,  2006; 
He,  2005;  He,  1998;  He,  1998;  He,  1999;  He,  1999;  He,  2000;  He,  2006;  He,  2005;  He,  2005;  He,  2005;  Ko<?ak 
et  al.,  2011;  Gepreel,  2011;  Cao  and  Bo.  Han,  2011;  Wazwaz,  2009).  Ezzati  and  Shakibi  (Ezzati,  2011)  solved 
the  Newell- Whitehead  equation  using  the  Adomian  decomposition  and  multi-quadric  quasi-interpolation 
methods.  They  concluded  that  the  Adomian  decomposition  and  multi-quadric  quasi-interpolation  methods  are 
reasonable  methods  to  solve  the  Newell- Whitehead  equation  with  acceptable  accuracy. 

In  the  present  work,  the  homotopy  perturbation  method  (HPM)  is  applied  to  obtain  the  closed  form 
solution  of  the  non-linear  Newell-Whitehead-Segel  equation.  Four  case  study  problems  of  non-linear  Newell- 
Whitehead-Segel  equations  are  solved  using  the  HPM.  The  trend  of  the  rapid  convergence  towards  the  exact 
solution  is  shown  when  compared  to  the  exact  solution.  The  Newell-Whitehead-Segel  equation  models  the 
interaction  of  the  effect  of  the  diffusion  term  with  the  nonlinear  effect  of  the  reaction  term.  The  Newell- 
Whitehead-Segel  equation  is  written  as: 
c?2ii 

— = k — - + au  - buq  (1) 

Ot  OX* 

Where  a?  b and  k are  real  numbers  with  k > and  q is  a positive  integer.  In  Eq.  (1)  the  first  term  on  the 

&u  y 

left  hand  side,  — , expresses  the  variations  of  with  time  at  a fixed  location,  the  first  term  on  the  right 

hand  side,  — - expresses  the  variations  of  with  spatial  variable  x at  a specific  time  and  the  remaining 

terms  on  the  right  hand  side,  ati  — , takes  into  account  the  effect  of  the  source  term.  In  Eq.  (1)  t)  is 

a function  of  the  spatial  variable  x and  the  temporal  variable  t,  with  X ■ ER  and  t >0.  The  function  u(:r,t) 

may  be  thought  of  as  the  (nonlinear)  distribution  of  temperature  in  an  infinitely  thin  and  long  rod  or  as  the  flow 
velocity  of  a fluid  in  an  infinitely  long  pipe  with  small  diameter.  The  Newell-Whitehead-Segel  equations  have 
wide  applicability  in  mechanical  and  chemical  engineering,  ecology,  biology  and  bio-engineering. 
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2.  The  Idea  Of  Homotopy  Perturbation  Method: 

The  homotopy  perturbation  method  (HPM)  is  originally  initiated  by  (He,  1999;  He,  2005;  Yildirim,  2009; 
He,  2006;  He,  2006;  He,  2006;  He,  2005;  He,  1998;  He,  1998;  He,  1999;  He,  1999;  He,  2000;  He,  2006;  He, 
2005;  He,  2005;  He,  2005;  Kogak  et  al.,  2011).  This  is  a combination  of  the  classical  perturbation  technique  and 
homotopy  techinique.  The  basic  idea  of  the  HPM  for  solving  nonlinear  differential  equations  is  as  follow; 
consider  the  following  differential  equation: 

= Oj-  (2) 

Where  E (it)  is  any  differential  operator.  We 

H(u,p)  = (1  -p)(F(tt)  - if (it))  + p(ff(u.))  = 0.  construct  a homotopy  as  follow: 


(3) 

Where  are  F(u),  H(u)  are  functional  operators  with  the  known  solution  v0.  It  is  clear  that  when  p is  equal  to 
zero  then  H(u,  0)  = F(u)  - H(u)  = 0,  and  when  p is  equal  to  1 , H(u,  1)  = E(u)  = 0.  It  is  worth  nothing  that  is  the 
embedding  parameter  p increases  monotonically  from  zero  to  unity  the  zero  order  solution  v0  continuously 
deforms  into  the  original  problem  E(u)  = 0.  The  embedding  parameter,  p £ [0,1],  is  considered  as  an  expending 
parameter(He,  1999).  In  the  homotopy  perturbation  method  the  embedding  parameter  p is  used  to  get  series 
expansion  for  solution  as: 

DO 

u = '^p,vi  = v0  + pv1+p2v2+p3v2-i (4) 

i =0 


When  p — M,  then  Eq.  (3)  becomes  the  approximate  solution  to  Eq.  (2)  as: 
u = va  + v1  + v2  + v3  H (5) 

The  series  Eq.  (5)  is  a convergent  series  and  the  rate  of  convergence  depends  on  the  nature  of  Eq.  (2)  (He, 
1999;  He,  2005;  Yildirim,  2009;  He,  2006;  He,  2006;  He,  2006;  He,  2005;  He,  1998;  He,  1998;  He,  1999;  He, 
1999;  He,  2000;  He,  2006;  He,  2005;  He,  2005;  He,  2005;  Kogak  et  al.,  2011).  It  is  also  assumed  that  Eq.  (3)  has 
a unique  solution  and  by  comparing  the  like  powers  of  p the  solution  of  various  orders  is  obtained.  These 
solutions  are  obtained  using  the  Maple  package. 

3 . The  Newell-Whitehead-Segel Equation : 

To  illustrate  the  capability  and  reliability  of  the  method,  four  cases  of  nonlinear  diffusion  equations  are 
presented. 


Case  I: 

In  Eq.  (1)  for  a = 2f  b = 3x  k = 1 and  q = 2 the  Newell-Whitehead-Segel  equation  is  written  as: 


at 


d1 


U 


— = — — - + 2u  - 3 u2 


dx2 

Subject  to  a constant  initial  condition 

u(x,0)  = X, 


(6) 

(7) 


We  construct  a homotopy  for  Eq.  (6)  in  the  following  form: 


— r-2u  + 3v2 


H(v,p)  = (1  - p)  — - 


~dv 

du0- 

+ p 

dv 

d2v 

dt  . 

dt 

dx 2 

= 0 


(3) 


The  solution  of  Eq.  (6)  can  be  written  as  a power  series  in  p as: 
v=va  +pv1  + p2v2  +■■■ 


(9) 
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Substituting  Eq.  (9)  and  Eq.  (7)  into  Eq.  (8)  and  equating  the  terms  with  identical  powers  ofp: 


dva 

dUfy 

£ 

O 

X 

O 

II 

pa: 

dt 

dt  ’ 

A, 

5 'ii. 

dun  d2vn 

P1 

± 

dt 

dt  dx2 

1^0  , 

£ 

H- 

X 

O 

II 

0, 

dv7 

d2v1 

+ u1(2  - 3u0)  , 

u2(x,0)  - 

p2: 

£ 

dt 

= dx 2 3"0"1 

0, 

d2v2 

— 3u2  +u2(2  — 3u0)  , 

p2: 

d 

1 ~dt 

= dx1  3“»U2 

^3  (t,0)  = 

0. 

Using 

; the  Maple  package  to  solve  recursive  sequences,  Eq. 

(10),  we  obtain  the  followings: 

= A, 

i.'j  (x,  t)  = A(2  — 3 A)  t, 

v2  (x,  t)  = 21(2  — 31)  (1  — 31)  ^ , (11) 

t3 

u3  ( x , t)  = 21(2  - 31)(2712  - 1S1  + 2)  — . 

By  setting  p = 1 in  Eq.  (9)  the  solution  of  Eq.  (6)  can  be  obtained  as 
u = u0  + v1  +■  v2  +■  +■  ...  Therefore  the  solution  of  Eq.  (6)  is  written  as: 

v (x ,t)  = 1 + 1(2  - 31) t + 21(2  - 31) (1-  31)  £ + 21(2  - 3l)(27lz  - 181  + 2)  ^ + 


The  Taylor  series  expansion  for 


--+A-Aeat 


is  written  as: 


-As 


2t 


■ I- 1 — As 


it 


= 1 + 1(2  - 31)  t + 21(2  - 31)  (1  - 31) 


2! 


+ 21(2  - 31)  (2712  -181  + 2)  — + ■■■ 

By  substituting  Eq.  (13)  into  Eq.  (12),  thus  Eq.  (12)  can  be  rewritten  as: 
-^le2t 

v (xrt  ) = - 

--+ A- Ae2t 


(12) 


(13) 


(14) 


This  is  the  exact  solution  of  the  problem,  Eq.(6).  Table  1 shows  the  trend  of  rapid  convergence  of  the 
results  of  50(^  t)  = (x,  t)  to  5 5 (xf  t)  = using  the  HPM.  The  rapid  convergence  of  the 

solution  toward  the  exact  solution,  the  maximum  relative  error  of  less  than  0.0046%,  is  achieved  as  shown  in 
table  1. 

Table  1 shows:  he  percentage  of  relative  errors  of  the  results  of  S0(r#  t)  = v0(r#t)  to 
55(x  0 = Sp=oVi  Of,  fc)  of  the  HPM  solution  of  Eq.  (6)  for  X=0.1 

Case  II:  In  Eq.  (1)  for  a — 1,  b =■  1,  k = 1 and  (7  — 2 the  Newell-Whitehead-Segel  equation  is 
written  as: 
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Table  1: 


percentage  of  relative  error  (%RE) 

t = 0.1 

Sa{x,t) 

0.1540788596 

Siix,  0 

0.01027226573 

S2(x,t) 

0.0002058041542 

Si  (x,t) 

0.00001949256084 

S4(xr  t) 

0.000001708760704 

S5{x,t) 

3.214500334  e-8 

t = .3 

S0(.x,t) 

0.3835101093 

S^x,  t) 

0.06910026503 

S2(x..t) 

0.003074197736 

S3(x,t) 

0.001358981068 

S4(_x,  t) 

0.0003091665977 

SB(x,t) 

0.00001230390524 

II 

Sa(x,t) 

0.4680703804 

S±(x,  0 

0.1063582391 

S2{x,t) 

0.005078839539 

S3(x,t) 

0.003988077959 

s4(x,  0 

0.001125247109 

Ss,(x,t) 

0.00004585924829 

du 

dt 


d2u 

= JP  + U~U 

Subject  to  initial  condition 


w(x,0)  = 


To  solve  Eq.  (15)  we  construct  a homotopy  in  the  following  form: 


H(vtp)  = (1-p) 


"3d 

du0- 

+ p 

dv 

d2v 

.dt 

dt  . 

dt 

dx2 

~ — 7 — u(l  - u) 


= 0. 


(15) 

(16) 


(17) 


The  solution  of  Eq.  (15)  can  be  written  as  a power  series  in  p as: 
v =va+pv1  + p2v2+--  (IS) 

Substituting  Eq.  (18)  and  Eq.  (16)  in  to  Eq.  (17)  and  equating  the  term  with  identical  powers  of  p , leads  to 

1 

u°(*'°)  = t TT7 

(l  + e^\ 


PQ-. 

d "2-Iq  5’W-  q 

dt  dt’ 

■H 

5 

dv1  duQ 

P ■ 

dt  dt 

dv2 

p . 

dt  dx2 

p3  = 

3li3  d2v2 

dt  dx2 

(19) 
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Using  the  Maple  package  to  solve  recursive  sequences,  Eq.  (19),  we  obtain  the  followings: 


uz(x,t)  = — 


(20) 


v3  (x,  t)  = 


By  setting  p — 1 in  Eq.  (18)  the  solution  of  Eq.  (15)  can  be  obtained  as 
v = + v±  + v2  + ^3  + — Therefore  the  solution  of  Eq.  (15)  is  written  as: 


/ 3/  -fVB 

f 1+8 % * j f 1+B'*  j 


The  Taylor  series  expansion  for 


s is  written  as 


(i  + .**)  (i  + .*)  3(i  + .*J 

125  (4(e^)2  - 


7e^  + 1 


216 


(l  + ^)B 


(21) 


(22) 


By  substituting  Eq.  (22)  into  Eq.  (21),  the  Eq.  (21)  can  be  reduced  to 
vQx,  t)  = 1 


(»M) 


(23) 


This  is  the  exact  solution  of  the  problem,  Eq.  (15).  Table  2 shows  the  trend  of  rapid  convergence  of  the 
results  of  5 0 (xr  t)  = Vq (:r , t)  to  = Zf=  q ^ t}  using  the  HPM  solution  toward  the  exact 
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solution.  The  maximum  relative  error  of  less  than  0.0008%  is  achieved  in  comparison  to  the  exact  solution  as 
shown  in  table  2. 


Table  2: 


percentage  of  relative  error  I 

(%RE) 

t=I 

II 

* 

x = l.S 

% = 1.8 

t = 0.1 

50(x,t) 

0.09374691950 

0.1010110875 

0.1051585828 

Site,  0 

0.003020748923 

0.003848096613 

0.004359470571 

s 2(x,t) 

0.00001102987604 

0.00002050070585 

0.00004357185889 

S3(x,t) 

0.000003008522247 

0.000003116144165 

0.000002949567767 

S4(x,  t) 

3.247853526  e-8 

7.597179054  e-8 

1.007266564  e-7 

55(x,t) 

2.475184340  e-9 

6.284562692  e-10 

4.103584300  e-10 

II 

U-j 

S0(x,t) 

0.2480807805 

0.2663112457 

0.2766522515 

5ife  0 

0.02225395729 

0.02841934500 

0.03220845585 

S2(x1t) 

0.0003852758255 

0.0003051199126 

0.0008095286291 

S 3(x,t) 

0.0002055820118 

0.0002152872468 

0.0002058262157 

S4(x,  t ) 

0.000005574634536 

0.00001431284116 

0.00001929380478 

S50t,t) 

0.000001472642781 

0.000001037290051 

5.725167943  e-7 

r-t 

II 

4^ 

0.3115489693 

0.3338572016 

0.3464690907 

s±(x,  t) 

0.03586206901 

0.04586954916 

0.05200200785 

S2(x,t) 

0.0009882331330 

0.0004901161293 

0.001569416719 

S3(x,t) 

0.0005982451688 

0.0006298759390 

0.0006050514464 

S4(x,  t) 

0.00001947940975 

0.00005317431566 

0.00007241822783 

Ss{x,t) 

0.000007711072685 

0.000005555853236 

0.000003218048969 

Table  2 shows:  the  percentage  of  relative  errors  of  the  results  of  50ti,  0 = i?0U,  fc)  to 
55(x  t)  = Gr,  t)  of  the  HPM  solution  of  Eq.  (15). 

Case  III:  In  Eq.  (1)  for  a = 1,  b = lfk  = 1 and  q = 4 the  Newell-Whitehead-Segel  equation 


becomes: 
du  d2u 
dt 


— = — ^ + u-  u 


dx 2 

Subject  to  initial  condition 


(24) 


BJf 

\1  + e-i™ 


2 


(25) 


We  construct  a homotopy  for  Eq.  (24)  in  the  following  form: 


= (i-p) 


~dv 

duQ- 

+ P 

dv 

ft. 

[dt 

dt  J 

dt 

1 

Cl, 

r-j 

1 

■V  + v4 


= 0. 


(26) 


The  solution  of  Eq.  (24)  can  be  written  as  a power  series  in  p as: 
v =vQ  + pv1  + p2v2  + ■■■ 


(27) 


Substituting  Eq.  (27)  and  Eq.  (25)  into  Eq.  (26)  and  equating  the  terms  with  identical  powers  of  p: 
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Qm 

^ dt  dt 


UQ(*,0)  = 


(l  + e-^y 


1 dv±  d2Ufy  du 


dt  dx 

dv2  e?2!^ 
dt  dx2 


3 _ 3u3  dzv2 


2 -^r+v bd-  vo3X 


p2^ 


dt 

^(J3^  + Vi_  , 


v1  (v,  0}  = 0., 

v2  ( x , 0)  = 0, 
u3  (x,  0)  = 0., 


(28) 


Using  the  Maple  package  to  solve  recursive  sequences,  Eq.  (28),  we  obtain  the  followings: 
a 

1 V 


'x'l  + 


vAx>0  = g 


3 

J — — i 


\ 


(l  + e--Uy 


/ 


uz  (x,  t)  = 


49 

50 


s \ 

-=i  1 


(t  + e^y  j 


(29) 


-27e^*+9^ 


v3  ( xr  t)  = 


343 

1000 


27e'JLCf  + 9 


(i+eU‘  / 


By  setting  p — 1 in  Eq.  (27)  the  solution  of  Eq.  (24)  can  be  obtained  as 
u = u0  + v1  + u2  + + ...  thus  the  solution  of  Eq.  (24)  can  be  written  as: 


vCr,t)=( ±3-)+^ 

\1  + eJH1)  5 


\ 


e-.'n> 


49 
t + — 

50 


+ 


343 

1000 


^26^*  - 3 j 


^i + y y y ^i + e ^ y j 

^4^6’^®*^  — 27e'^x  + 9^e,>fi°*^  ^_3 

V (i+^)T  j 


(30) 
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2 

The  Taylor  series  expansion  for  ^ tank  ^ ^ J + 


is  written  as: 


Comparing  Eq.  (31)  with  Eq.  (30),  thus  Eq.  (30)  can  be  reduced  to 


2, 


(32) 


This  is  the  exact  solution  of  the  problem,  Eq.  (24).  Table  3 shows  the  trend  of  rapid  convergence  of  the 
results  of  50(^t)  = to  5s(^t)  = using  the  HPM  solution  toward  the  exact 

solution.  The  maximum  relative  error  of  less  than  0.007%  is  achieved  in  comparison  to  the  exact  solution  as 
shown  in  table  3.  We  can  conclude  that  the  HPM  is  one  the  most  suitable  and  friendly  method  in  solving  the 
Newell-Whitehead-Segel  equation. 

Table  3 shows:  the  percentage  of  relative  errors  of  the  results  of  5.j(x  t)  to 

t)  = Ef=0i^  Gx  t)  of  the  HPM  solution  of  Eq.  (24). 


Case  IV: 

In  this  case  we  will  examine  the  Newell-Whitehead-Segel  equation  for  a = 3>  b = 4,  k = 1,  q = 3, 
chi  d2u 
dt 


— = — — ^ + 3 u - 4 u , 


dx 2 

Subject  to  initial  condition 


(33) 


We  construct  a homotopy  for  Eq.  (33)  in  the  following  form: 

H(vrp)  = (1-p) 

The  solution  of  Eq.  (33)  can  be  written  as  a power  series  in  p as: 
Li  = V±  + p2V2  H 


~dv 

dua- 

+ p 

dv 

d2v 

3t 

dt  J 

dt " 

dx2 

= 0. 


(34) 

(35) 

(36) 


Substituting  Eq.  (36)  and  Eq.  (34)  into  Eq.  (35)  and  equating  the  terms  with  identical  powers  of  p : 
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,dv1=du1 

p ■ st  dt' 


A _ 


P 


dVl  d2v0  dvQ 

■it 


(r,  0)  = 0,, 


2 c?u2  3zu1 

1 2u  4-  Tii 

u2< 

12d2  D(,  d3( 

:x,0) 

x,0) 

= o„ 
= 0. 

^ dt  dx 2 

3 . 9u3  _ ^2 

LjL'Uq  eJ-^  | J1  L'^ 

-1-  ^11  1 ?7I 

^ dt  dx 2 

Table  3: 

T J L1-?  JL^cJqcJ^ 

percentage  of  relative  error  (%RE) 

X = 1 

x = l.S 

x = 1.8 

t = 0.1  50Ut) 

0.09323514238 

0.1045372618 

0.1099789256 

Sxix,  t) 

0.001725301473 

0.003516790607 

0.004498961488 

S2{x,t) 

0.0002100221036 

0.0001216772460 

0.00005170881962 

S3(x,t) 

0.00001173612982 

0.00001539927439 

0.00001472293353 

S4(x,  t) 

3.620778386  e-7 

2.929797495  e-7 

5.746720664  e-7 

s 5 0,0 

5.414213656  e-8 

3.614791446  e-8 

1.434500497  e-8 

t = .3  50Otr,  t) 

0.2396054415 

0.2702774140 

0.2852064054 

St(_x,  t) 

0.009390569965 

0.02330928345 

0.03106724554 

S2(x>t) 

0.005215736776 

0.003375969333 

0.001825357483 

S3(x,t) 

0.0007262150389 

0.001037577469 

0.001023346798 

S4ix,  t) 

0.00009555513293 

0.00004044412369 

0.0001029640746 

S5(x,t) 

0.00003280557979 

0.00002473085989 

0.00001198739229 

II 

a- 

0.2962707966 

0.3351845099 

0.3542406943 

Si(x,  t) 

0.01219208796 

0.03518330953 

0.04811460219 

S2(x,t) 

0.01183961035 

0.008037424274 

0.004713588643 

S3(x,t) 

0.001990820701 

0.002987593293 

0.002996129616 

S4(xr  t) 

0.0004128329726 

0.0001164720814 

0.0003682018367 

55(x,0 

0.0001681121669 

0.0001337455766 

0.00006941923143 

Using  the  Maple  package  to  solve  recursive  sequences,  Eq.  (37),  we  obtain: 


u0O,t)=  - — 

-Jsx 


-J&x 


VifoO  =2 


(eV6je  + e = *J 


t. 


, , 81  3 

uz(x,0  = — T 


e^eV 


+e»*)tz 


( 1±\ 

/eV6*  + g j 
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v3  (x,  t)  = 


243  3 


a ' 


16  J 4 


|eV&ae  + e a 


By  setting  p — 1 in  Eq.  (36)  the  solution  of  Eq.  (33)  can  be  obtained  as 
v = vQ  + v±  + v2  + v3  + ... 

Thus  the  solution  of  Eq.  (33)  can  be  written  as, 

v ft,  t)  = 


i?  . > I?  , , 3i 

^ .T'+(.^)’  t(.f J 


243 

16 


(39) 


The  Taylor  series  expansion  for 


1 r-  (fx-lt) 


is  written  as 


V fcJt’ 


_|_  g V 2,  * 3 7 


Vfoe 


3 e 

4 _ _ 

gV&s  -[-  g *a  * 


9 

+ T 


_ V t> 


+ 


31  +.¥“),» 
4 j4  . _ 

3 V&Je 


243 

TfT. 


^e^e~ 
3 

4_ 


_|_  e : 

T*  + (e^)2  + ^eT*j  ^ 


-|-  g 1 

By  substituting  Eq.  (40)  into  Eq.  (39),  the  Eq.  (39)  can  be  reduced  to: 


(40) 


v(xf  t)  = 


,V6*  -|_g(s*  /) 


(41) 


This  is  the  exact  solution  of  the  problem,  Eq.  (33).  Table  4 shows  the  trend  of  rapid  convergence  of  the 
results  of  50(^t)  = v^{x},t)  to  Ss{xft)  = using  the  HPM  solution  toward  the  exact 
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solution.  The  maximum  relative  error  of  less  than  0.037%  is  achieved  in  comparison  to  the  exact  solution  as 
shown  in  table  4. 


Table  4: 


percentage  of  relative  error  (%RE) 

x = .2 

II 

x = .8 

t = 0.1 

s0O,t) 

0.1591055626 

0.1376714465 

0.09890297594 

Si(x,  t) 

0.007038530707 

0.009754764500 

0.01176933513 

S2(x,t) 

0.002482893305 

0.001788283111 

0.0004607861238 

S3(x,t) 

0.0001958736131 

0.0002690488417 

0.0002512962487 

s4(.x,  t) 

0.00004554598541 

0.00002344226507 

0.00001229642845 

S5{x,t) 

0.000004998678372 

0.000006399733729 

0.000003681928211 

t = . is 

Sa(.x,t) 

0.2155132753 

0.1864801193 

0.1339670591 

S±(_x,  0 

0.01698529379 

0.02214248709 

0.02558157620 

S2(x!t) 

0.007422698254 

0.005232453593 

0.001127443849 

S3(_x,t) 

0.001011675135 

0.001318032849 

0.001182316216 

S4(.x,  t ) 

0.0003016921205 

0.0001450263132 

0.0001001951821 

■S' 5 ft,  0 

0.00005638430243 

0.00006875985734 

0.00003732434558 

t = .2 

■S' ,3  Or,  0 

0.2605557486 

0.2254546357 

0.1619662983 

■S'lOf,  0 

0.03164334146 

0.03938244479 

0.04388750886 

s2(x,t) 

0.01561926168 

0.01076040185 

0.001819030418 

S3{x,t) 

0.003225422844 

0.004022798070 

0.003478947590 

s4  Of,  0 

0.001110362592 

0.0004931283715 

0.0004433744920 

s'bOvO 

0.0003119309814 

0.0003646043520 

0.0001870028440 

Table 

4 shows: 

the  percentage 

of  relative  errors  of  the 

results  of  t)  = 

55(x  t)  = Gr,  t)  of  the  HPM  solution  of  Eq.  (33). 


4.  Conclusion: 

In  the  present  work  the  exact  solution  of  the  Newell- Whitehead- Segel  nonlinear  diffusion  equation  is 
obtained  using  the  HPM.  The  validity  and  effectiveness  of  the  HPM  is  shown  by  solving  four  non-homogenous 
non-linear  differential  equations  and  the  very  rapid  convergence  to  the  exact  solutions  is  demonstrated 
numerically.  The  trend  of  rapid  and  monotonic  convergence  of  the  solution  toward  the  exact  solution  is  clearly 
shown  by  obtaining  the  relative  error  in  compared  to  the  exact  solution.  The  rapid  convergence  towards  the 
exact  solutions  of  HPM  indicates  that,  using  the  HPM  to  solve  the  non-linear  differential  equations,  a reasonable 
less  amount  of  computational  work  with  acceptable  accuracy  may  be  sufficient.  Moreover  it  can  be  concluded 
that  the  HPM  is  a very  powerful  and  efficient  technique  which  can  construct  the  exact  solution  of  nonlinear 
differential  equations. 
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